Spin-entropy driven melting of the charge order in Lai 5 Sr .5CoO4. 
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We studied the melting of the charge order in the half-doped tetragonal perovskite Lai.sSi'o.sCoCU 
by elastic neutron scattering. We found that diffuse peaks, corresponding to the breathing-type 
modulation of oxygen positions in the C0O6 octahedra, disappear above T c = 825(27) K. This melt- 
ing of the diffuse superstructure is reversible (no change in the correlation lengths upon annealing 
is observed) and accompanied by a large nonlinear thermal expansion along the tetragonal c-axis, 
with a cusp at T c . We conclude that quenched disorder of the doped Sr ions is at the origin of this 
anisotropic charge glass state. We also suggest that its melting is driven by the transition from the 
intermediate- to the high-spin state of Co 3+ ion, which is favored by the spin entropy at high T. 



PACS numbers: 71.28.+d 71.45.Lr 75.40. Gb, 

Understanding the cooperative charge ordered phases 
in doped transition-metal oxides is a key problem in 
the physics of colossal-magnetoresistance (GMR) mate- 
rials and high-temperature superconductors Eva Electron 
correlations in these systems are driven by contribu- 
tions from several equally important interactions, such 
as Hund's coupling, splitting of the on-site electronic 
levels by the crystal field, Jahn- Teller (JT) distortion, 
Coulomb repulsion between the conduction electrons, su- 
perexchage, double exchange, etc. Competing ground 
states lead to a rich variety of phases, depending on the 
electronic configuration of the 3d metal ion, its crystal en- 
vironment, the doping level, and so on. To understand, 
and ultimately predict, the nature and the physical prop- 
erties of a given transition-metal oxide at a particular 
doping, it is extremely important to study the origin of, 
and the relation between, charge order (CO), spin order 
(SO) and orbital order (00) in different compounds. 

The case of half-doping, which is relevant for the CMR 
materials, is of special interest for several reasons. First, 
the spatial structure of the CO instability is particu- 
larly simple - it is a checkerboard-type arrangement of 
the doped charges in the tetragonal plane. Second, it 
is argued to be of most technological relevance in per- 
ovskite manganites.B Finally, an orbital order, concomi- 
tant with CO, was recently directly observed in sever 
half-doped manganites by X-ray resonant scattering! 
Consequently, a number of theoretical papers appeared, 
aimed at explaining the origin of the observed CO/OO at 
half-doping, its relation with-snin order and implications 
for the transport properties na To verify these theories, 
and eventually achieve clear understanding of this impor- 
tant problem, more experimental input is vital. 

In a recent neutron scattering studya we character- 
ized spin and charge order in the half-doped cobaltate 
Lai.5Sr .5CoO4. We found that the spins order at T s w 
30 K in a quasi-two-dimensional incommensurate glassy 
structure, while checkerboard short-range CO persists to 
at least ~ 600 K and shows no measurable anomaly 



at T s , i.e. the freezing of charge and spin correlations 
are independent instabilities, inherent to the electronic 
and crystal structure of Lai.5Sro.sCo04. Both magnetic 
and structural scattering were similar to those observed 
in the isostructural manganite Lao.sSri.sMnO^EHI where 
CO and SO transition temperatures are much closer, 
and corresponding instabilities are often presumed to be 
closely related. The main difference in our case is the 
much shorter coherence range: correlations within the 
tetragonal a-b planes extend only over £ s - 80 A for SO, 
and less than 30 A for CO, while £ co > 200 A in 

Lao.5Sr1.5MnO4.Ej Despite its short-range glassy nature, 
CO in Lai.5Sro.5Co04 has no less drastic consequences 
for transport properties than that in manganitejij re- 
sulting in an insulating state, and an activation behavior 
of the electrical conductivity with E a ~ 6000 K.E3 In 
fact, such an apparently large activation energy (i.e. a 
weak temperature dependence of the conductivity), is ob- 
served in the pseudo-cubic perovskite cobaltites RC0O3- 
R=La,Pr,Nd,Sm,Eu,Gd at temperatures 400 - 600 K,E3 
where a semiconductor-to-metal cross-over occurs upon 
heating. Comparative studies of the isostructural Co and 
Mn compounds give a unique opportunity to understand 
the relative importance of the anisotropic crystal field 
energy, which in the cobaltites can be larger than the on- 
site Hund's coupling, and the inter-site spia-Lnteraction, 
which is often quenched for the Co 3+ ions.ljO 

Here we focus our attention on the nature of the 
CO glassy state and the charge melting process in 
Lai.5Sr .5CoO4. We study the temperature evolution of 
the CO correlation lengths and the charge order param- 
eter up to and above T c , and investigate the effect of an- 
nealing and slow cooling of the sample through the CO 
transition. Such procedure could distinguish between the 
cooling-rate-dependent disorder of the mobile ions (non- 
stoichiometric oxygen) or charge carriers- which can re- 
order into an equilibrium arrangement ,Ld and that re- 
lated to the quenched disorder of the immobile Sr 2+ ions. 

We studied the same high quality single crystal sam- 
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pie of Lai.5Sro.5Co04 as in Ref. At all tempera- 
tures the crystal remained in the tetragonal "HTT" phase 
(space group IA/mmrri), with low T lattice parameters 
a = 3.83 A and c = 12.5 A, Fig. 1(a). In Ref. | we used 
the "orthorhombic" indexing, based on the space group 
FA/mmm with a unit cell that is two times larger. This 
accounts for breaking of the translational symmetry in 
the a-b plane by the checkerboard CO superstructure, 
and is convenient for describing the spin order. Heija, 
however, we will use the "HTT" IA/mram notation,E2l 
which is relevant for the charge-disordered phase, and is 
more appropriate for characterizing the CO instability. 
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FIG. 1. (a) Crystal structure of Lai.sSro.sCoO^, arrows 
show the directions of oxygen displacements £o(i),o(2) accom- 
panying the CO modulation. (b),(c) Filling of the electronic 
levels in (b) Co 2+ with S=3/2, and (c) in Co 3+ in the LS 
(S=0), IS(S=1) and HS (S=2) states. 

Experiments were performed on BT2, BT4 and BT9 
thermal beam 3-axis neutron spectrometers at NIST 
Center for Neutron Research. PG(002) reflections were 
used at the monochromator and analyzer, supplemented 
by two PC filters to suppress the contamination from 
the higher order reflections. The energy of the scattered 
neutrons was fixed at Ef = 14.7 meV, beam collima- 
tions were 80' - 42' - 62' - 100' on BT2 and BT4, and 
40' - 56' - 51' - 100' on BT9. The cylindrical sample was 
mounted either in a displex refrigerator or in a 850 K 
vacuum furnace (BT9) with the axis vertical, allowing 
wavevector transfers in the (hhl) reciprocal lattice plane. 
Some scans on BT9 were also done at ambient in the 
(hkO) plane. Sample mosaic is <■ 0.25° in the (hhl) plane, 
and < 0.4° in the (hkO) plane. Normalization of the 
scattering intensity was performed using the incoherent 
scattering from a vanadium standard. 

An extensive survey of elastic scattering from the 
super-structure induced by the CO is presented in Fig. 
2. It consists of broad commensurate peaks centered at 
Q = (h + 0.5, k + 0.5, 1) reciprocal lattice units (rlu) with 
h, k, I integer, which corresponds to a checkerboard ar- 
rangement of the Co 2+ /Co 3+ valence in the a-b plane. 
An increase in the intensity at large Q, characteristic of 
the scattering from a structural modulation, is apparent. 
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FIG. 2. Elastic scattering from the diffuse structural modu- 
lation accompanying charge order in Lai.sSro.sCoOj. (a)-(f) 
scans collected on BT2 at T = 6 K in the (hhl) reciprocal 
plane. (g),(h) scans collected on BT9 at ambient T « 295 
K in the (hkO) reciprocal plane. Horizontal bars show the 
full width at half maximum (FWHM) of the spectrome- 
ter elastic resolution along the scan direction. Solid curves 
through the data are the resolution corrected fits to the sin- 
gle peak cross-section (^). Shaded peaks show the net inten- 
sity which would result from four incommensurate peaks at 
Q c = {(0.5 ± e, 0.5 ± e, /), (0.5 ± e, 0.5 =p e, I)} with e = 0.016. 

In all scans the widths of the peaks are substantially 
larger than the instrument elastic resolution, which re- 
veals the short-range nature of the CO. The peaks 
are broader in the c-axis direction, so that the super- 
structure is an anisotropic glass with the correlation ex- 
tending by about 4 lattice repeats in the a-b plane, and 
only 2/3 of a lattice repeat along the c axis, i.e. only be- 
tween the two nearest C0O4 planes. The intensity mod- 
ulation with / shown in the Fig. 2 (d), (e), exhibits a 
rather peculiar long-periodic pattern, which translates 
into a modulation of a small object in real space. In 
Ref. |^ we showed that scattering from a breathing-type 
distortion of the oxygen octahedra surrounding the Co 
ions provided a good account of a smaller set of data. 
This involves displacements of the in-plane oxygens e^'L , 
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and apical oxygens Eq( 2 )' a l° n g the corresponding Co-0 
bonds, as shown in Fig. 1 (a). The neutron scatter- 
ing elastic cross-section arising from such a modulation, 
short-range-periodic in the crystal lattice, has the form 
of the factorized "lattice Lorentzians" ,113 



%{q) = NI DW \F{q)\ 2 i- £ Ct (q) (q) Cf (q) , (1) 

sinh^" 1 



cosh 1 — cos (q ± Q) ■ a 



etc. (2) 



Here F(q) — J2^{q ■ e tl )b fl e~ 1 (Q' r ^ , N is the number 
of unit cells in the crystal, Idw is an effective Debye- 
Waller intensity prefactor, is the displacement of the 
atom /i inside a unit cell from its nominal position r^, 
is its scattering length, Q is the modulation wavevector, 
£a = Cfc = £ft ano - = are the in-plane and inter-plane 
correlation lengths measured in lattice units (lu), and n± 
is the number of terms in the last summation. 

We performed a global fit of all measured intensities 
at each T with Eq. ([!]), corrected for the instrument 
resolution, and implying an equal-weight superposition 
of the commensurate modulations with Q = (0.5, 0.5, Z) 
for I = 0, f . This corresponds to an equal probability 
for the in-phase and anti-phase stacking of the checker- 
board CO structure for neighboring planes along c. At 
low temperatures the Debye- Waller factor was fixed at 
Idw = 1 ; an d the displacements of the O nuclei from 
the nominal positions ro(\) = (0.5, 0, 0), (0, 0.5, 0) lu and 
r o(2) = (0,0,0.173) lu, the correlation lengths £h and 
the background from the q- independent incoherent 
elastic scattering and the Bragg scattering by the poly- 
crystalline aluminum in the sample environment were 
refined. Adding a modulation of the La/Sr position, 
rLa.Sr = (0.5,0.5,0.141) lu, as well as varying Zo(2) and 
ZLa.Sr, did not improve the quality of the fit. The solid 
lines through the data in the panels of Fig. 2 result from 
Eq. (|l]) with parameters obtained in the global fit and 
summarized in the Table 1. Clearly, the scattering from 
the oxygen displacements gives a very good description 
of all measured intensities. 

In the Ref. [| we also found magnetic elastic scattering 
which appears below T rj 30 K at slightly incommen- 
surate positions Q m = with e w 0.016 
and I odd-integer. If this incommensurability of the 
spin order is the result of a stripe order of the doped 
charges, related to a small deviation of the La:Sr ra- 
tio from the nominal 1.5:0.5, or some non-stoichiometric 
oxygen, we would also expect four incommensurate CO 
peaks at Q c = {(0.5 ± e, 0.5 ± e, ± e, 0.5 T e, I)}, 
as observed in the doped cupratesO However, for the 
scans in the (h,h,l) plane, Fig. 2 (a)-(f), the major 
contribution of the out-of-plane peaks, accepted by the 
coarse vertical resolution of the spectrometer and cen- 
tered at h = 0.5, would hinder the experimental obser- 
vation of the incommensurate splitting.0 To clarify this 
issue we made a number of scans through the CO scat- 



tering around (1.5,1.5,0) and (2.5,0.5,0) in the (h, k, 0) 
plane. No splitting was apparent in any of these scans 
either, as illustrated by two representative datasets in 
Fig. 2 (g),(h). The overall peak width is always domi- 
nated by the short correlation length, which smears the 
implied small splitting. In other words, the correlation 
range of the CO, £h « 4 lu, is so small compared to the 
wavelength of the expected incommensurate modulation 
A e = 1/e w 60 lu, that the latter is undetectable. The 
cross-section resulting from four incommensurate peaks 
with e = 0.016 is shown by the shaded areas in the Fig. 
2. It gives an almost identical fit to the measured inten- 
sities, with only a slightly larger ^ to "compensate" for 
the splitting, see Table 1. 
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FIG. 3. Change in the diffuse super-structural scattering 
concomitant with melting of the CO in Lai.sSro.sCoO,^ 

Figure 3 shows the evolution of the diffuse CO scat- 
tering upon heating. Although a small peak may still 
be present in the /i-scan at T — 809 K, we find that 
it completely disappears at 850 K. Each measurement 
at high T was preceded by a sample re-alignment, since 
the melting of the charge order in Lai.5Sr .5CoO4 is ac- 
companied by an anomalously large thermal expansion 
along the c-axis, Fig. 4(a). We find that the temper- 
ature dependence c(T) saturates somewhere above 800 
K, with a kink indicative of a second order phase tran- 
sition. This supports our suggestionEI that the decrease 
in the c lattice spacing below T c is proportional to the 
charge order parameter r\ c . Fitting it with the expression 
c(T) = c(850 K) - a(T c - Tf we get T c = 825(27) K and 
j3 = 0.59(15), indistinguishable from the mean-field value 
0.5. The in-plane lattice constant a, on the other hand, 
stays almost unchanged, with some tendency to decrease 
upon heating, and exhibits no anomaly at T c . These find- 
ings clearly indicate that electrons involved in the CO are 
localized in the out-of-plane d 3z 2_ r 2 orbitals. 

Interestingly, to within our errors the width of the peak 
in h does not change upon heating, i.e. the in-plane cor- 
relation length stays constant, and CO melting mostly 
involves the loss of the inter-plane coherence with increas- 
ing T, Fig.4 (b). In the fits of the scans for T > 300 K we 
fixed the oxygen displacements and varied the "Debye- 
Waller" prefactor to account for the decrease in the CO 
scattering intensity, which follows the disappearance of 
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FIG. 4. Temperature dependence of the structural param- 
eters related to the CO in Lai.sSro.sCoO^. (a) change of the 
lattice constants relative to the low-T values, (b) in-plane 
(left scale) and inter-plane (right scale) correlation length, for 
T > 650 K they were fixed at £ h = 3.5 lu, £j = 0.32 lu. (c) 
intensity ( "Debye- Waller" ) prefactor. Broken curves are the 
guides for the eye, solid lines are the critical-law fits (see text). 

the order parameter rj c , Fig. 4(c). The solid line in Fig. 
4(c) is the best fit to the expression Idw(T) ~ (T c — T) 7 , 
which gives T c = 821(48) K, identical with that refined 
from the c(T), and 7 = 0.92(11). The exponent 7 agrees 
well with 7 = 2/3, expected from low ~ Vc' an( i a ^ so 
with the mean-field value 7=1. 

We found that upon slow cooling, even after annealing 
the sample for a few hours at T = 850 K, the CO peaks 
reappear below T c with the same width and intensity as 
before. The reversible nature of the melting transition is 
important for understanding the origin of the charge glass 
phase in Lai.sSro.sCoC^. It suggests that the quenched 
disorder of the dopant Sr ions, which strongly interact 
with the doped charges, is the most probable cause for 
the loss of the CO in-plane coherence, and not a non- 
stoichiometric oxygen or rapid quench through T c . 

A key for understanding the Co 3+ /Co 2+ mixed valence 
phases and cobaltites in general is provided by the pecu- 
liar structure of the 3d electron levels. In a cubic crys- 
tal field 3d levels are split into a lower ti g triplet and 
an e g doublet of d x 2_ y 2 and d 3z 2_ r 2 orbitals. In cobalt 
the Hund splitting of spin "up" and "down" 3d levels 
is anomalously small, and though Co 2+ (3d 7 ) is in the 
"normal" S=3/2 state t\ g e 2 g , Fig. 1(b), Co 3+ (3d 6 ) may 
be found in either the high t\ g e 2 g (HS, S=2), intermediate 
t\ g e\ (IS, S=l), or low spin state t% (LS, S=0), Fig. 1(c). 
If LS were the ground state, as in LaCoOs, a decrease in 
the free energy F = E — Tln(2S + 1) due to the higher 
paramagnetic entropy may drive transitions to IS and 



HS states with increasing temperatureJialia On the other 
hand, it is clear from Figure 1(c) that only the IS S=l 
state is JT-active, i.e. favors tetragonal splitting of the 
e g doublet, as it has one electron there. In either the HS 
or LS state, such splitting gains no energy because both 
e g levels are equally occupied. 

We conclude, that Co 3+ ions in Lai.sSro.5Co04 are in 
the IS state at low T, which favors the JT-distorted CO 
phase. At higher T, a spin-entropy driven transition to 
the HS state occurs (similar IS-HS transition is observed 
in LaCo03 at « 600 Kt§) with consequent disappearance 
of the JT modulation and conspicuous melting of the 
charge order. Therefore, it appears that local effects such 
as the spin-state transition and JT level splitting are very 
important and quite efficient in driving the CO phase in 
Lai.5Sro.5Co04, and, perhaps, in other materials. 
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